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Abstract 

We study total muon capture rates in light (A ~ 6-18) nuclei, taking into 
account renormalizations of the nuclear vector and axial vector strengths. We 
estimate the influence in the results of uncertainties of the spin-isospin interac- 
tion parameter g' and nuclear densities. A few of these reactions are theoretical 
benchmarks for physics involving searches for neutrino oscillations. New exper- 
iments in muon capture in several targets are suggested, in the light of some 
discrepancies with theory, crudeness of some experimental results and relevance 
to neutrino physics. 
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Since its discovery in the forties [0, nuclear muon capture (NMC) has been a 
probe par excellence for both weak interaction studies and exploring nuclear dynam- 
ics @. In the context of QCD, the nucleus provides a new vacuum, wherein quark and 
gluon condensates can have entirely different values from those in the free hadrons 
@, thereby opening the possibility of nuclear renormalization of the hadronic weak 
couplings UJ]. On the other hand, particle-hole and/or isobar- hole excitations provide 
a more traditional theoretical setting for describing the nuclear renormalization dy- 
namics ||. We explore here the latter in NMC via the study of inclusive (or total) 
capture rates. 

The importance of NMC in light nuclei (here, those with mass number A < 20) 
transcends nuclear physics. Thus, muon capture reactions in 12 C and 16 0, to take two 
examples, provide us with opportunities to test the nuclear dynamical models, useful 
in interpreting the neutrino reactions, for example, those studied by the KARMEN 
H and the LSND [|7| collaborations. The latter have bearings on the question of the 
neutrino oscillations, thus, physics beyond the standard model. The inclusive NMC 
rates in these targets are important benchmarks for the theory of nuclear response, 
of great value in setting the scale of our understanding within the standard model. 
While inclusive NMC rates are generally difficult theoretically to compute || ^ [10], 
11 , HJ, these are experimentally relatively easy to measure ||12|| , particularly in muon 



factories, where high muon purity and a decent stopping rate are readily available 
T3fl . Regrettably, such experiments have been few and far between in the recent 



times ||12|| . Our theoretical work here hopes to elicit new experimental interest in this 
subject, of importance to nuclear physics and particle physics. 
The reaction of our interest is 

fi-(lS) + (A,Z)^u, + (A,Z-iy, (1) 

wherein the nucleus (A, Z—l)*, formed in the reaction, is not experimentally detected. 
The rate A c of the inclusive reaction is experimentally inferred from the disappearance 



rate of the muon B 12 



Theoretically, three common ways of calculating A c are: (1) by computing and 
summing exclusive channels [pT0| , (2) by using the sum rules || and (3) by using the 



elementary particle method [II]. We follow a different method |J, which assumes the 
reaction (1) proceeding as 

IA-(lS) + \p]^u„ + [n], (2) 

where [p] and [n] are nucleons belonging to two different local Fermi seas, characterized 
by the Fermi momenta kp^(r) (i = p,n), given by the medium densities Pi(r) = 
/c^j(r)/37r 2 . Then, the differential rate for the reaction (2) is given, in the local 
density approximation, by |§ 

^~EEI^%?1, (3) 

where we suppress in Eq. (3) phase space factors and measures of integration, shown 
fully in Ref. [9]; |T| 2 is the square modulus of the transition amplitude, suitably 
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renormalized in the nuclear medium, summed and averaged over the initial and 
final states respectively; the function U{qo,q) is the nuclear Lindhard function for 
the particle-hole excitation. The three-momentum transfer q is fixed by the outgoing 
neutrino momentum. In our approach the energy transfer, go, to the nucleus is given 

by 

q = E fl + Q' -Q- E„. (4) 
Here is the available energy from the muon 

£ M = rrv - El s , (5a) 

where the second term on the right-hand side of (5a) is the muon binding energy. Q' 
comes into play in the nucleus with N ^ Z, 

Q' = E nF - E PF ; (56) 



Q is the threshold for the nuclear reaction (1) to begin 



Q = M f - Mi , (5c) 

the nuclear mass difference for the target and daughter nuclei. For light nuclei under 
consideration in this Letter (target nuclei between G Li and 18 O), the Q value varies 

from a low of 0.67 MeV , in the case 

of 14 N fL, Uq transitions, to a high of 21.15 MeV, 
in the case of the 14 C target, metastable, but still experimentally accessible, because 

of its very long life, leading to 14 C 14 B transitions; indeed, Q values exceeding 
10 MeV are quite common in light nuclei. We find a significant effect of the Q value 
on the total capture rate in light nuclei, since the right-hand side of (3) is a sensitive 
function of the available energy in the reaction, (2). 

The calculations in [|10j use a RPA approach to nuclear structure. Kolbe et al. 
use a continuum RPA which sums up over the excited nuclear states above nucleon 
emission. Auerbach et al. use a standard RPA, including pairing, that allows one to 
include the contribution of all final states. In both cases the RPA evaluation reduces 
the results with respect to the calculation with the single particle orbitals, but the 
reduction is stronger in the calculations of Auerbach et al. 

Our approach might look simplified with respect to the ones just mentioned, but in 
fact it is also an RPA approach built up from single particle states of an uncorrelated 
local Fermi sea. This method in practice is a very accurate tool when the excitation 
energy is sufficiently large such that relatively many excited states contribute to the 
process, and in particular if a large fraction of it comes from excitation to the contin- 
uum, as it is the case in capture. The adaptation of the method to finite nuclei via 
the local density approximation has proved to be a rather precise technique to deal 
with inclusive photonuclear reactions [T5| , response functions in electron scattering 



16| and deep inelastic scattering fl7| . 



Obviously, because of its nature, the method only applies to inclusive processes, 
summing over relatively many final states and it is not meant to evaluate transitions 
to discrete states. 
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The simplicity of the method, however, allows for improvements over the tradi- 
tional methods: 1) the sum SS|T 2 | is done relativistically; 2) the nucleon momentum 
of the Fermi sea is considered in the calculations which usually is taken equal to 
zero; 3) the Z e // approximation is avoided and instead our evaluation uses directly 
the muon wave functions. Nonrelativistic muon wave functions are used, but the 
relativistic effects were evaluated in || and amounted to about 2%. The fi~ binding 
energy is also considered. 4) In addition, our approach builds up the RPA correlations 
by allowing the Ah excitation on top of the ph excitation considered in [TD|. This 



is a major difference with the other theoretical RPA approaches mentioned Jl0|. A 
Landau- Migdal interaction is used to account for the propagation of these ph and Ah 
states in the nucleus. 

The Ah excitation leads to extra quenching of the transition strength with respect 
to the one found by Kolbe et al. The role of this Ah excitation in the quenching of 
Gamow Teller transitions was emphasized by Rho and later by Brown and Rho in ref 
This quenching could be interpreted as an effective quenching of gA, since this 
coupling term gives the largest contribution to the /i capture rate; we can, instead, 
talk in terms of quenching of the transition strength. The pseudoscalar current is also 
renormalized, but this term plays a smaller role in /i capture. 

Our method was used in ||, where a good overall agreement with experiments 
over the periodic table was found. The method has proved even more accurate than 
anticipated in the ref. |J, when the improvements of the present paper are done. The 
basic framework of our approach has been described in detail in ||. We stress here 
only significant improvements of this approach: (1) We use a new Lindhard function 
|T8|] , which takes into account a shift of the neutron and proton Fermi seas due to 



the experimental Q value |T4[]. Its imaginary part coincides with the the results of 
the Eq.(3) in terms of the ordinary Lindhard function, but its real part is changed 
slightly with respect to the old one. The new Lindhard function avoids pathologies 
of the ordinary Lindhard function in the limit {qo,q) — > 0. This limit, however, 
does not occur in the muon capture kinematics, but it does play an important role 
in the neutrino scattering at low energies. (2) The experimental Q value and the 
theoretical one, Q', are used in the evaluation as indicated in eqs. (3), (4). The 
Q-value in general is important in determining the NMC rate in light nuclei. (3) The 
uncertainties coming from the Landau-Migdal parameter g' are included by varying 
g' in the range [[HJ 

(j 7 = 0.7 ±0.1. (6) 

This is an important parameter influencing the nuclear response in our approach. 
This parameter occurs in the particle- hole and the Delta-hole interactions @, ||. (4) 
In many nuclei, the radial form of the nuclear density has parametric uncertainties. 
We take them also into account in our calculation by letting the density parameters 



vary within the experimental errors of |20 



In Table I, we display our calculated total capture rate A c , obtained by integrating 
(3), for g' = 0.7 and for one set of standard nuclear density parameters pCfl . The 
present approach has significantly improved the agreement between the theory and 
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experiment, compared with results reported in ||. 

The principal uncertainty of our theoretical estimate of the NMC rate stems from 
that of the parameter g' . For the range of g', given in equation (6), this results in an 
uncertainty of about ±(5 — 7)% in the NMC rate for all the nuclei considered in this 
paper and reported in table I. Thus, to take a physical example, that of NMC by 12 C, 
we vary g' from 0.6 to 0.8 and this results in a variation of the rate from 3.78 x 10 4 s~ l 
to 3.43 x 10 4 s _1 , a ±5% effect around the central value of 3.6 x 10 4 s -1 , corresponding 
to g' = 0.7. Our experience with many electroweak and strong processes points 



towards a value of g' = 0.7 [15] , |T6| , F7| . Similarly, other implicit parameters of the 
particle-hole or the Delta-hole interactions are determined from other processes. 

The nuclear densities p PjTl , enter in our calculation via the local density approx- 
imation |J. There are three cases, for A equal to 12, 14 and 16, in which radial 
uncertainties are readily investigated in our approach by calculating the capture rate 
for many sets of density parameters determined from electron scattering fl20| . They 



produce relatively small uncertainties of ±(2 — 3)%. The theoretical capture rates 
quoted for these nuclei in table 1 should be read with this additional uncertainty due 
to the nuclear densities. Thus, in 12 C, the total capture rate varies from 3.45 x 10 4 s _1 
to 3.60 x 10 4 s _1 , keeping the spin-isospin parameter fixed. This is a ±2 % variation 
due to uncertainty of the nuclear density around the central value. Assuming the 
uncertainties of g' and nuclear densities to be independent, we have a theoretical un- 
certainty in 12 C NMC rate of ~ 6%: this is a conservative estimate of our theoretical 
error. Thus, the l2 C NMC rate is, according to our theoretical calculation, 

A C ( 12 C) = (3.60 ± 0.22) x 10 4 s" 1 . (7) 



This is in good agreement with the best experimental determination so far 
taking the world average of the best determinations, with their errors in quadrature: 

A e c xp { 12 C) = (3.80 ± 0.10) x 10 4 s- 1 . (8) 

In the following, we compare our work with the recent work of Kolbe et al. and 
Auerbach et al. for this reaction. 



Kolbe et al. [K| have done a continuum RPA calculation, which can compute 
with a reasonable accuracy the NMC rates to the particle-unbound states in l2 B. 
Using their calculation of the 1994 paper [10], we obtain the inclusive capture rates 
for the 12 C nucleus, wherein the particle bound states are not included for the excited 
nucleus (see their Table I). We get theoretical rates for two different potentials to be 
3.42 and 3.34 in units of 10 4 s _1 . These numbers are slightly above the experimental 
value for the capture to the particle unbound states quoted by Kolbe et al. Taking a 
conservative lower limit of 7000 for the muon capture rate to the particle-bound 
states directly from the experiment, we can translate the calculation of Kolbe et al. 
as 

A C ( 12 C) ~ (4.04 - 4.12) x 10 4 s _1 . (9) 
Auerbach et al. have several variations of their HF-RPA model and they obtain 
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A C ( 12 C) ~ (3.09 - 3.64) x 10 4 s _1 . (10) 

Recent researches in neutrino physics have brought NMC in 12 C in sharper fo- 
cus as a theoretical benchmark. One has to keep in mind the results of the recent 
(u e , e~) and experiments || [/]], along with the results of some recent theoret- 



ical calculations p0|, [21|| . The latter reaction is the inverse of the NMC. In particular, 
the recent LSND experiment has produced an inclusive (i/^, yr) rate, which is in dis- 



agreement with the calculation of Kolbe et al. [10| by a factor of about 1.5. However, 



Auerbach et al. ||10|, in a recent paper find that this disagreement can be eliminated 



in their Hartree-Fock- Random-Phase Approximation (HF-RPA) calculation that in- 
cludes pairing. With the present approach, we obtain a rate which is within 15 % of 
the experimental result when the maximum uncertainties in the theoretical as well as 
experimental results are taken into account pl| . 

In the neutrino physics context, the NMC produces the following situation: while 
Kolbe et al. obtain a result for NMC rate which is within 5% of the experimental 
result and explain the (y e ,e~) inclusive cross sections rather well, they overestimate 
the inclusive (z/^/i - ) cross section by (45-50) %. In the calculations of Auerbach et 
al. it is possible in one version of their model to obtain a NMC rate (A c = 

3.64 x 10 4 s _1 ) in fair agreement with the experiment, but this version of the model 
predicts slightly higher values for the (z/ M , /i~) inclusive cross sections and a large 
value for the {y e , e~) inclusive cross sections when compared with the experiments . 
On the other hand, the version of the model, which can explain the inclusive neutrino 
reactions , predicts a NMC rate of 3.09 x 10 4 s _1 which is smaller than the experimental 
value. Our present method, when applied to l2 C reproduces the experimental results 
of muon capture as well as those of the (z/ e , e~) inclusive cross sections quite well, but 
it also overestimates the experimental values for the (z/ M ,yU~) inclusive cross sections 
by about (15-20)% This underscores the important role of the NMC in selected 
targets as benchmark reactions to calibrate the nuclear theory, which is, in turn, used 
to delineate standard model physics. From this point of view, the NMC rates in 12 ' 13 (7 
and 16 are particularly important. 

The above discussion, comparing results of different accurate methods, gives us 
an idea of present theoretical uncertainties in calculations of weak nuclear reactions. 

We now come to a very subtle physics in NMC, the isotope effect. In Table II, 
we display our calculations for the ratios of the rates of NMC for a pair of isotopes. 
We compare here our results with experimental ones wherever available, and the 
prediction of the well-known formula of Primakoff ||: 

where 5 ~ 3.15. The Primakoff formula yields Pauli blocking of the muon capture 
process rather approximately. 

From Table II, we can see that our calculation does much better than the Primakoff 
formula in three out of four cases where there are data on isotope shifts. For example, 
our calculation yields for 13 C/ 12 C and 18 0/ m O ratios 0.82 and 0.78, in fair agreement 
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with the experimental values of 0.90 ± 0.01 and 0.86 ± 0.15 respectively, in contrast 
to the predictions of the Primakoff formula of 0.71 and 0.59 in the two cases. Our 
approach is able to take into account the special properties of individual target nuclei 
such as the shell closure in 16 and the large Q-value needed in 14 C. No such physics 
is included in the relatively simple formula (11). 

Detailed calculations, such as those in JIDJ, should likewise do better than the Pri- 
makoff formula. They also take into account the role of the experimental Q values. 
The case of 7 Li / e Li does not come up too good in our approach. This is under- 
standable, since, in very light nuclei, the local density approximation can become less 
accurate. 

We now make some remarks on the experimental situation and future prospects 
at the muon factories. It is important to have muon capture rates in light nuclei 
measured with high precision. In several targets, old || and newer [T^] experiments 
disagree quite strongly (e.g. 9 Be, 14 N, ie O etc.), even though the more recent exper- 
iments are of higher statistical accuracy. While our theoretical estimates agree with 
most of the recent determinations of the A c , there are significant disagreements in 
7 Li and 14 iV. Some targets, 6 Li and 7 Li, have relatively poorly known A c and de- 
serve better experimental determination. 14 C, being a radioactive nucleus, will pose 
experimental challenges to prepare a dense enough nuclear target. However, it is very 
interesting theoretically: the Q-value for NMC in this target is the largest in the 
nuclei studied here, with an enormous effect on the rate. Hopefully, the long half-life 
of 14 C would make an experiment with it feasible. 

Our worst agreement with experiment comes in 7 Li and 14 N targets. We predict 
the inclusive NMC rates 3.4 x 10 3 and 8.7 x 10 4 s _1 in these two cases, while the most 
accurate experiments give (2.26 ±0.12) x 10 3 s _1 and (6.93 ±0.008) x 10 4 s _1 respec- 
tively [T2|]. Given the quality of our agreement in other targets with experimental 



results, this disagreement invites an experimental reconfirmation, with a special care 
on the systematic errors of the experiments. Further experimental studies are also 
needed to investigate possible non-statisticality in the hyperfine states of the muonic 
atom , in particular, in the case of 14 N, due to hyperfine conversion. 

In summary, we have studied total NMC rates in the framework of a theory 
wherein vector and axial-vector strengths in nuclei are appreciably renormalized due 
to particle-hole and Delta-hole correlations. Large effects from the consideration of 
the Q- values of the NMC are seen, in some nuclei. A reasonable (± 0.1) variation 
of the Landau-Migdal parameter g' around its preferred value 0.7 translates into a 
±(6 — 7)% uncertainty of the NMC rate for the nuclei studied in this paper. Similarly 
the nuclear radial uncertainties in some of these nuclei contribute at most a ±(2 — 3)% 
variation of the capture rate. 

Since renormalizations of weak vector and axial-vector strengths, studied here, 
also occur in the electromagnetic and strong processes, we should interpret them as 
those for nuclear strengths rather than for weak nucleon couplings in nuclei. This 
has important consequences in our understanding of the nuclear QCD effects, in 
particular, the Bjorken sum rule in nuclei. 

Precise theoretical understanding of the NMC in 12 C, 13 C and w O would provide 
benchmarks in low-energy neutrino physics, with important bearings on the issue of 
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neutrino oscillations and physics beyond the standard model. Although small, there 
are still defferences between the predictions of different accurate models and it would 
be most desirable to examine and test these differences in precise experiments in 
future. 
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Table I: Total muon capture rates in s -1 in various target nuclei, compared with 
the most accurate experimental results cited in the literature, wherever available [|i~2|| ; 
A dash in the last column means that no data are available in that case; in case of 
two or more experimental results of comparable quality, we display the value with 
best precision here. We use g' — 0.7 and one standard set of radial densities (20 



Target 


This calculation 


Experiment 


G Li 


4.68 x 10 3 


(4.68 ±0.12) x 10 3 


7 Lt 


3.41 x 10 3 


(2.26 ±0.12) x 10 3 


9 Be 


8.84 x 10 3 


(7.4 ±0.5) x 10 3 


io B 


2.67 x 10 4 


(2.78 ±0.07) x 10 4 


U B 


1.86 x 10 4 


(2.19 ±0.07) x 10 4 




3.60 x 10 4 


(3.76 ±0.04) x 10 4 


13 C 


2.95 x 10 4 


(3.38 ±0.04) x 10 4 


14 c 


2.34 x 10 4 




U N 


8.67 x 10 4 


(6.93 ±0.008) x 10 4 


is at 


6.34 x 10 4 




16 


1.16 x 10 5 


(1.026 ±0.006) x 10 5 


17 


1.06 x 10 5 




1SQ 


9.00 x 10 4 


(8.80 ±0.15) x 10 4 



Table II: Isotope effect (with g' = 0.7) calculated by us, compared with the 



Primakoff formula (Eq. (14)) and experiment |L2 . 



S 



Isotope pair 


This calculation 


Primakoff 


Exp. 


7 Li/ 6 Li 


0.73 


0.47 


0.48 ± 0.03 


n B/ 10 B 


0.79 


0.67 


0.79 ± 0.03 


n C/ l2 C 


0.82 


0.71 


0.90 ± 0.01 


1 A - — -j / 1 O - — x/ 

14 cy 12 c 


0.65 


0.47 


- 


14 cy 13 c 


0.79 


0.66 


- 


15 jV/ 14 iV 


0.73 


0.75 




17 0/ 16 


0.91 


0.78 




180/160 


0.78 


0.59 


0.86 ± 0.15 


180/170 


0.85 


0.75 
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